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We present an efficient implementation of the Density Matrix Renormalization Group (DMRG) 
algorithm that includes an optimal ordering of the proton and neutron orbitals and an efficient 
expansion of the active space utilizing various concepts of quantum information theory. We first 
show how this new DMRG methodology could solve a previous 400 KeV discrepancy in the ground 
state energy of ®®Ni. We then report the first DMRG results in the pf + g9/2 shell model space for 
the ground 0^ and first 2"*" states of ®^Ge which are benchmarked with reference data obtained from 
Monte Carlo shell model. The corresponding correlation structure among the proton and neutron 
orbitals is determined in terms of the two-orbital mutual information. Based on such correlation 
graphs we propose several further algorithmic improvement possibilities that can be utilized in a 
new generation of tensor network based algorithms. 


Large scale Shell Model calculations in physically 
sound valence spaces are the prime choice in the nu¬ 
clear spectroscopy of light and medium mass nuclei. The 
catch is that the dimensions of the basis (the number 
of Slater determinants) grow as the product of the two 
combinatorial numbers made with the number of single 
particle states (n,l,j,m) in the valence space and the 
number of particles, for neutrons and protons. Modern 
shell model codes can cope with dimensions ©(fO^^jP. 
Going beyond this limit requires clever truncation and 
extrapolation schemes. In addition, different approxima¬ 
tion methods may be implemented. Among them, per¬ 
haps the most accurate is the Monte Carlo Shell Model 
(MCSM) that stochastically samples the Hilbert space 
to find the relevant basis states for the description of a 
particular eigenstat^. Several improvements have been 
added to the original methodology along the years mak¬ 
ing the MCSM a robust and accurate technique to study 
medium-mass nuclei. The most important are sequen¬ 
tial conjugate gradient to select the Slater determinants 
and the parity and angular momentum projection of 
the determinants, plus an energy variance extrapolation 
methocP. 

Another efficient numerical tool to approximate the 
exact wave function in a truncated basis, is the DMRG 
method^ that was earlier introduced in nucl ear structure 
either in the particle-hole basis (phDMRGjP^ or in the 
j-coupling scheme (JDMRG)^. Both methods found dif¬ 
ficulties to treat systems beyond the limits of an exact 
diagonalization. However, phDMRG in the j-coupling 
scheme has been successfully applied to the Gamow shell 
model of weakly bound lighPH nuclei due to the weak 
entanglement between the valence space and resonances 
with the discretized continuum. Moreover, a standard 
implementation of DMRG to nuclear structure found a 


serious discrepancy of 400 KeV for the ®®Ni ground 
state energy in the pf shelP^. These results hindered 
further application of the DMRG method to medium 
size nuclei. Meanwhile, there have been important ef¬ 
forts to implement DMRG in quantum chemistry (QG- 
DMRGjpi^by utilizing various concepts of quantum infor¬ 
mation theorjli^ti^. Nowadays, DMRG is capable to pro¬ 
vide the low-lying energy spectrum of complex molecules 
with great accuracj®"2o| and it has ranked among the 
standard multireference QC methods. 

In this letter we will make use of the new DMRG tech¬ 
niques developed in QC to overcome first the 400 KeV 
discrepancy found in the ®®Ni calculations in the pf shell. 
Afterwards, we will show how DMRG can deal with nu¬ 
clei well beyond the limits of an exact diagonalization by 
studying ®^Ge in the pf and in the enlarged pf + g9/2 
valence space, and compare our results with benchmark 
calculations of MCSM. As a byproduct of the method, 
we will depict the single-site entropy and mutual infor¬ 
mation of ®'^Ge, shedding new light into the landscape of 
entanglement and correlations in nuclear structure. 

In our DMRG implementation we study the most gen¬ 
eral Hamiltonian with one- and two-body interaction 
terms given as 

LI — ^ ^ 2 ^ ^ L)^/37(5CqC^C,5C^ , (1) 

ot OLj3^5 

where c\ and Cq creates and annihilates a particle with 
quantum numbers a = (n, l,j, m, Tz)- 

In the so-called representation an orbital can be 

either empty or occupied, thus the dimension of the local 
Hilbert space, Ai of a single orbital, q = dim is 2. The 
full Hilbert space of a hnite system comprising N orbitals, 
AW, is built from tensor product spaces of local orbital 
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spaces Ai, which can be written as = 0^iAi. 

In the DMRG method, quantum correlations are taken 
into account by an iterative procedure that variationally 
minimizes the energy of the Hamiltonian given by Eq. Q . 
The method eventually converges to the full Configu¬ 
ration Interaction (Cl) solution within the selecte d ac- 
tive orbital space. In the two-site DMRC varianlP^, 
aw is approximated by a tensor product space of four 
tensor spaces defined on an ordered orbital chain, i.e., 
■^DM^RG “ ® ^i +1 ® Aij ^2 ® The basis states 

of the comprises i orbitals to the left of the chain 
(Z = left ) and comprises N — i — 2 orbitals to the 
right of the chain (r = right). These states are deter¬ 
mined through a series of unitary transformation based 
on the singular value decomposition (SVD) theorem by 
going through the ordered orbital space from left to right 
and then sweeping back and forttP^*211, xhe number of 
block states, Mi = dim and Mr = dim re¬ 

quired to achieve sufficient convergence can be regarded 
as a f unctio n of the level of entanglement among the 
orbit alJi^IlIl xhe maximum number of block states 
Mmax = max {Ml, Mr) required to reach an a priory 
defined accuracy threshold, is inherently determined by 
truncation error, dcTR, when the Dynamic Block State 
Selection (BBSS) approach is usecP^. During the initial 
sweeps of the DMRC algorithm the accuracy is also af¬ 
fected by the environmental error, dsgweep- The latter 
error can be reduced significantly by taking advantage of 
the Cl based D ynami cally Extended Active Space proce¬ 
dure (CI-DEAS)P^^and using a large number of DMRC 
sweeps until the energy change between two sweeps is 
negligible. In the CI-DEAS procedure the active space 
of orbitals is extended dynamically based on the orbital 
entropy prohl^Sl_ M^ax depends strongly on the orbital 
ordering along the on e-dimensional chain topology of the 
DMRC methocPMSl. xhere exist various extrapolation 
schemes to determine the truncation-free solutiorP^. In 
this work, due to the high level of entanglement of the 
nuclear wave functions, we carry out the extrapolations 
as a function of the total number of block states M. Us¬ 
ing E{M) = E{M ^ oo)+ XiM^^ we have estimated 
the truncation-free solution where E{M —>■ oo),Ai and 
X 2 are free parameters of our fit!^. We have performed 
between 30 to 90 sweeps by requiring the energy change 
between two sweeps to be below 10“^ MeV. 

The amount of contribution to the total correlation 
energy of an orbital can be quantified by the single-orbital 
von Neumann entropjff^, Si = —Tipi Inp^ where pi is the 
reduced density matrix of orbital i. The two-orbital von 
Neumann entropy sij is constructed similarly using the 
reduced density matrix, pij of a subsystem built from 
orbitals i and j, and the mutual information = Sy — 
Si — Sj describes how orbitals are correlated wi th each 
other as they are embedded in the whole system.E^E^IIIl 
The orbital ordering is determined by the minimization of 
the entangle ment distance expressed as /dist = lij 
where rj > l. b^d4| -y^rork we used 77 = 2 in order to 

carry out the optimization task using concepts of spectral 


graph theorjBH. 



FIG. 1. (Color online) The DMRG ground state energy for 
®®Ni and ®^Ge in the p/-shell and ®^Ge in the pf+g9/2 valence 
space, are shown as a function of 1/M. Diamonds (DMRG), 
Crosses (Exact Diagonalization), Circle (MCSM). The solid 
lines are the fits that produce our extrapolated value, while 
the dotted lines joining the points not used in the fit are just 
to guide the eye. 

Our DMRC code for nuclear structure calculations is 
composed of two phases. The preprocessing phase in 
which the orbital ordering and the active space vectors 
are optimized by calculating the one-orbital entropy and 
two-orbital mutual information using a hxed small num¬ 
ber of block states. The production phase in which an 
accurate calculation is performed using large fixed num¬ 
ber of block states or the DBSS procedure in order to 
reach an a priory set accuracy threshold The pre¬ 
processing phase takes only a small fraction of the total 
computational time. 

First we perform the DMRC calculation for the sd nu¬ 
cleus ^®Si employing the USD interactiorP^. This system 
was treated almost a decade ago by the DMRC methocP^. 
By keeping fixed M = 1024 block states we have re¬ 
produced the exact value of the 0+ ground state energy 
and excitation energy of the first 2+ state. Using the 
quantum information entropy based ordering optimiza¬ 
tion, the CI-DEAS procedure, and the DBSS approach 
with Mmin = 256 and fcxR = 10“^ — 10“® we deter¬ 
mined the ground state energy within the a priory set 
accuracy threshold. This drastic improvement clearly 
demonstrates the importance of the quantum entropy 
based optimization procedures. 

Next w e con sider the nucleus ®®Ni with the KBS 
interactioiP^m. The exact results in M valence space 
were obtained by the code ANTOINhW The basis di¬ 
mension in m-scheme for Jz=0 is 1.1 x 10®. In Fig. [^a) 
the DMRC ground state energy for ®®Ni is shown as a 
function of 1/M. For small M values a downward cur- 
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vature governs the behavior of the scaling function (in¬ 
dicated by the dotted lines) while for large enough M 
values, an inflection point is reached and an upward cur¬ 
vature becomes apparent. Our lowest variational energy 
obtained with M = 4096 states is E^g^M = 4096) = 
—78.451 MeV, thus the error compared to the exact 
value, ilexact = -78.465 MeV, is 1.4 x lO-^ MeV. 

Using the fit function defined above we have esti¬ 
mated the ground state energy in the M —> oo limit as 
E^g\M —)■ oo) = —78.463(4) MeV. The deviation from 
the exact ground state energy is in the order of a few 
keV’s. Therefore, using entropy based optimized DMRG 
methodology we could solve the previous 400 keV dis¬ 
crepancy in the ground state energy of ^®Ni reported in 
Ref. do]. 

Next we study ^"‘Ge in the pf valen ce sp ace for which 
we have used the GXPFIA interactiorPS®. The DMRG 
results for the pf valence space are displayed in Fig.[^b). 
Our lowest variational energy obtained with M = 4096 
states is E^\m = 4096) = —304.163 MeV which ex¬ 
trapolates to E^\m —>• oo) = —304.182(3) MeV while 
the exact energy obtained with the code ANTOINE is 
A'exact = —304.183 MeV. Therefore, the error is again in 
the keV range. 

Being confident with the results obtained in the pf 
shell, we proceed to explore the ability of DMRG to 
describe accurately the structure of nuclei in valence 
spaces that exceed the limits of an exact diagonaliza- 
tion. With this in mind, we study ^'^Ge in the extended 
space pf+g9/2 (dimension in m-scheme 1.7x 10^"* ^) wh ich 
was already considered using the MCSM methocpESl. In 
this case the GXPFIA interaction has to be completed 
with the matrix elements between the p/-shell orbits and 
the g9/2 obtained in a standard G-matrix calculation. 
In addition, the spurious center of mass contamination 
is treated with the Lawson’s ansatP^l^SI Within this 
extended valence space we found a significantly slower 
convergence rate as a function of M as can be seen in 
Fig. s c). Even using up to M = 3072 block states the 
inflection could not be reached, thus the upward curva¬ 
ture did not become apparent yet. Therefore, the esti¬ 
mated energy using our fit function overshoots the ref¬ 
erence MGSM energy (Amcsm = —306.066) MeV and 
can provide only a lower bound which we found to be 
E^\m —)• oo) = —307.4 MeV. An upper bound can 
be estimated using a second order polynomial fit giv¬ 
ing E^\m —>■ oo) = —305.5 MeV. The extrapolated 
MCSM reference energy lies within the two bounds given 
above. Due to the slow scaling of the energy as a func¬ 
tion of 1 /M significantly more block states are needed to 
provide a reliable extrapolation. 

Besides optimization procedures the entanglement 
analysis is also very important to obt ain physic al in¬ 
formation encoded in the wave function^l^l^IMISI!, xhe 
single-orbital entropy and two-orbital mutual informa¬ 
tion obtained with DMRG for ®^Ge are shown in Fig. 

The single orbital entropy profiles reflect the strong 
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FIG. 2. (Color online) Single-site entropy and mutual infor¬ 
mation obtained with DMRG for ®'*Ge in the pf -|- g9/2-shell. 
The mutual information matrix elements are shown on a lad¬ 
der topology with time-reversed pairs in the rungs. Circles 
and diamonds label proton and neutron orbitals respectively, 
and sites are denoted by l,j,m with +m outside the ladder 
and —m inside. 


configuration mixing in the wave function. The mag¬ 
netic states belonging to the same orbital have the same 
entropy (within the numerical accuracy). The entropies 
of the pl/2 and p3/2 orbitals are very close to the upper 
limit ln2 = 0.693 which corresponds to the maximally 
mixed state. Orbitals with largest entropies contribute 
the most to the correlation energy thus their accurate 
treatment is mandatory. As expected, the entropy pro- 
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files show that the pl/2, p3/2 and /5/2 orbits play a 
dominant role in the physics of the system, whereas the 
presence of the /7/2 orbit is non-negligible and that of 
the 59/2 is minor. Notice that these three highly en¬ 
tangles orbits close a pseudo-SU(3) symmetry which en¬ 
hances the quadrupole-quadrupole correlations as it will 
be apparent in the mutual information diagram. 

The two-orbital mutual information shows how orbitals 
are correlated with each other. The rungs of the ladder 
display the degree of entanglement between time reversed 
states of like particles. The fact that the mutual informa¬ 
tion is approximately equal for the pl/2, p3/2 and /5/2 
orbits, and independent of their projections, is con¬ 
sistent with the presence of a strong T=1 proton-proton 
and neutron-neutron pairing coherence. The same co¬ 
herence, though significantly less intense, is seen in the 
/7/2 and the g9/2. Proton-neutron T=1 pairing correla¬ 
tions could also be seen between time reversed and charge 
conjugated states for the pl/2, p3/2 and /5/2 orbits. We 
could also see signihcant entanglement between proton- 
neutron maximally aligned states for the p3/2 and /5/2 
orbits, which could be related to J = 2jz pairing and/or 
quadrupole-quadrupole in the T=0 channel. Similarly, 
quadrupole correlations can be observed inside the lad¬ 
der (proton-proton and neutron-neutron) as well as those 
connecting opposite sites (proton-neutron) of the ladder 
for the pl/2, p3/2 and /5/2 orbits. 

Let’s consider now the excitation energies of the 2+ 
states. It is worth to remark that even in the case when 
the ground state energy could not be determined with 
the desired accuracy, the energy difference between two 
states can be obtained more accurately due to cancela¬ 
tion of errors by targeting the two lowest lying eigenstates 
simultaneously. The exact value of the 2“'" excitation en¬ 
ergy for ^®Si can be reproduced with M = 1024 block 
states. For ®®Ni, using fixed M = 2048 block states, we 
obtain AE = 5.218 MeV to be compared with the ex¬ 
act value 5.125 MeV, with an absolute error of 93 keV. 
For ^"‘Ge in the p/-shell we get AE = 0.922 MeV com¬ 
pared with the exact value 0.906 MeV, the absolute er¬ 
ror being now 16 keV. Using the fitting procedure to 
extrapolate the 2+ excitation energies to M=oo we ob¬ 
tain 5.121(3) MeV and 0.907(2) MeV for ^SNi and ^^Ge. 
Indeed the agreement is excellent. 

When the p/ -I- g9/2 valence space is considered we 
found a slower convergence rate as a function of 1/M 
due to the more complex entanglement structure in the 
system. Fig. [^a) and (b) show the 2+ excitation ener¬ 
gies of ®^Ge as a function of 1/M for the p/-shell and 
p/ + 59/2-shell, respectively. The solid lines represent 
the fitted values. In the latter case the exact solution is 
not available, thus the extrapolated MCSM energie^^, 
AE — 0.919 (without reordering) or AE = 0.890 (with 
reordering) can be taken as a benchmark reference and 
are indicated by the blue diamond symbol. Our extrap¬ 
olated excitation energy AE = 0.90(2) is in good agree¬ 
ment with the MGSM result. 

In this work we have demonstrated that the DMRG 



FIG. 3. (Color online) The energy gap between the 0 and 
2 + states for ^’^Ge with pf and pf + g9/2 valence space are 
shown as a function of 1/M. Circles (DMRG), Cross (Exact 
Diagonalization), Diamond (MCSM). The solid lines are our 
fits. 


method including various novel optimization algorithms 
based on quantum information theory can be applied 
efficiently to nuclear structure calculations for medium- 
mass nuclei in extended valence spaces with dimensions 
exceeding the limits of an exact diagonalization. The 
DMRG results for the ground state energy of ®'^Ge in the 
pf + g9/2 valence space could be further improved in or¬ 
der to obtain a reliable estimate for the truncation-free 
limit. However, the excitation energy of the first 2+ is 
in excellent agreement with MCSM benchmark results. 
By calculating the single-orbital and two-orbital mutual 
information we have determined the correlation struc¬ 
ture among the orbitals. These entanglement graphs, 
based on recent developments in QC applications^, con¬ 
stitute a novel and highly precise tool to picture the wave 
function correlation properties. The analysis of these 
entanglement graphs suggests new alternatives for effi¬ 
cient truncation methods that can be developed on gen¬ 
eral grounds. A straightforward extension of the DMRG 
algorithm could be the use of local tensors with Ai of di¬ 
mension 5 = 4, but unlike the QC implementations, the 
states to include in this tensor should be time-reversed 
pairs in order to optimize the treatment of nuclear pair¬ 
ing correlations. As a final remark, the use of tree-tensor 
network state (TTNS), which is a recent development in 
quantum information theorjI^lHini could tackle the prob¬ 
lem of having equal single entropy values for each orbital 
group that makes inefficient a sequential treatment of 
these states within an ordered chain. Within this scheme, 
several orbitals with equally large entropy values could 
form the central shells of the TTNS network. These new 
developments could certainly open the field of medium to 
heavy nuclei to highly precise spectroscopic calculations. 
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